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The sensitivity of experiments looking for annual modulated signals is discussed and analyzed. The 
choice of energy intervals for rate integration and enhancing the signal-to-noise ratio of the pre- 
dicted WIMP signal is addressed. Special emphasis is placed on quantifying the minimum required 
exposure, MT, under experimental conditions. The difficulty reduces to establish the proper sepa- 
ration between the rate due to the unmodulated part of the WIMP signal and the rate of spurious 
background present in any experiment. The problem is solved by placing an upper bound to the 
unmodulated part of the signal using the best exclusion plot. We find that the lowest backgrounds 
achieved result in exposures in the range MT ~ 5 — 100 kg yr for masses > 100 GeV depending 
on the energy threshold of the detector. While our results are valid for Ge and Nal detectors, the 
formulae derived apply to other elements as well. Prescriptions are given to estimate exposures for 
higher background experiments. 
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I. INTRODUCTION 

Weakly interacting massive particles (WIMP's) are 
most likely the main constituents of the galactic halo. 
Constraints on the strength of the interaction and mass 
of the WIMP's have been obtained by several experi- 
ments using germanium [|, silicon [|, and Nal(Tl) |] 
crystals. For a coupling constant of the order of the weak 
coupling constant, the mass of these WIMP's must be, in 
any case, heavier than the 40 GeV lower bound already 
imposed by the measurement of the width of the . 
The same measurement has also placed bounds on some 
supersymmetric candidates. 

Drukier et al. suggested that a distinctive signa- 
ture of WIMP's would be the annual modulation of the 
detection rate. The modulation originates in the orbital 
motion of the Earth around the Sun that produces a vari- 
ation in the relative velocity between the Earth and the 
WIMP's, thus altering the dark matter flux during the 
year. Extracting this fluctuation from data that, cer- 
tainly, contain the contribution of spurious events is a 
task that demands care. The problem to be confronted 
with is the following: to extract a feeble signal from 
the data, we require a large signal-to-noise ratio so that 
the relative statistical fluctuations are small compared 
with the sought after signal. At the same time, efforts 
to reduce the background of unwanted (no dark mat- 
ter) events prevent achieving large signal-to-noise ratios. 
There is no compromise solution possible; if a low count- 
ing rate has been attained, a long exposure is required to 
detect an even smaller modulation. 

The main purpose of this work is to examine this is- 
sue quantitatively and to establish some guidelines upon 
which upcoming data, like those from the Sierra Grande 
experiment j^, can be analyzed. Although the subject 
of annual modulation, and its details, have been dis- 
cussed by several authors we feel that the sensi- 



tivity of modulation searches, namely the amount of ex- 
posure needed to resolve a given fluctuation embedded in 
a background, has not been yet addressed quantitatively 
in the literature. A semi-quantitative approach was ad- 
vanced in Ref. |^2| where results for required exposures 
were presented. Our procedure to obtain minimum ex- 
posures however, relies on simpler assumptions and give 
lower estimates. 

In the first section we briefly introduce the formalism 
needed for calculating the expected signal-to-noise ratios. 
For the sake of simplicity, we illustrate our calculations 
for the case of a germanium detector later on we ex- 
tend the results to a sodium iodide detector. We take 
into account the influence of the crossing of the spectra 
due to their annual fluctuation and choose energy inter- 
vals of integration that maximize the expected signal- 
to-noise ratios of the modulated rate. The next section 
analyzes quite generally the sensitivity of this type of 
experiments, addressing its difficulties and common mis- 
uses. Estimates of the exposures for the best experiment 
are presented together with guidelines for experiments 
with higher backgrounds. Most of our results are not 
strongly influenced by astrophysical uncertainties. 

II. PREDICTED SIGNAL 

The expected total rate of events due to the recoil of 
nuclei elastically scattered by WIMP's will be the prod- 
uct of the cross section, the WIMP's flux, and the number 
of target nuclei in a detector of atomic mass number A. 
For an incoming WIMP of mass and velocity v the 
differential counting rate in the recoil-energy interval T, 
T + dT , is given by 

dR Na ( Phalo\ \ '^^ ( \^ in ^\ 
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where Na is Avogadro's number, phaiojnix num- 
ber density of WIMP's, and /(u, CTi,, Vg) is the Maxweh- 
Bohzmann velocity distribution of the WIMP's for an 
observer on the Earth The velocity distribution is a 
function of Vg (t) , the velocity of the Earth with respect 
to the galactic rest frame varying annually and ay, the 
dispersion velocity of the WI MP's in the galactic halo. 
The integration limits of ( |2.l| ) are, Vmax, the maximum 
velocity of the WIMP's {vmax — Vs + Vesc, where Vesc 
is the escape velocity from our galaxy), and Umm, the 
minimum WIMP velocity necessary to contribute to a 
particular energy of the recoil spectrum. 

The differential cross section in Eq. (^]^) is given, in 
general, by 
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dT 
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where ao{fTi^) is the cross section at zero momentum 
transfer (of a heavy Dirac neutrino, for example), Tmax = 
iji^v^ /rriM (/i is the WIMP-nucleus reduced mass and 
toat the nucleus mass), and F'^{T) is a standard Bessel 
nuclear form factor |T^] that depends on the momentum 
transfer and which takes into account the loss of coher- 
ence of the interaction. The calculated recoil spectrum 
must also be convoluted with a relative efficiency func- 
tion for germanium p^ ] accounting for the efficiency of 
the recoiling process in generating an ionization signal. 
We refer therefore to the deposited energy E, instead of 
T, in the following formulae. Table | lists the values of 
Phaioi Vesc, o'v, and Vq , the velocity of the Sun around 
the galactic center, used in this work. 



The figure explores also the dependence of the crossing 
energy on different halo parameters. The halo density 
for example, has no effect on the crossing energies since 
it affects only the absolute rates. The crossing is neither 
affected by the assumed escape velocity, Wesc, or the value 
of Vq (which in turn modifies Vg) but, it is sensitive to the 
dispersion velocity, a^, which for WIMP masses higher 
than 500 GeV changes by 2 keV, when a^ is increased 
from 270 km to 300 km s~^. The dispersion velocity 
of WIMP's in the galactic halo is a rather uncertain pa- 
rameter, ranging from 246 km s~^ to 323 km s^^ [Q, so 
here we adopt a standard value of 270 km s^^. 
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FIG. 1. Crossing energies in Ge as a function of WIMP 
mass for two different dispersion velocities. The dotted line 
corresponds to a dispersion velocity (j„ = 300 km s~^ and the 
dashed line to (j„ = 270 km s^^. 



Crossing energies 

It is known that, the differential energy spectra of 
June and December cross at specific energies that de- 
pend, among other factors, on the WIMP mass The 
reason for the crossing of the spectra is that in June the 
velocity distribution increases by 15 km s~^ enabling the 
WIMP's to deposit a higher fraction of their kinetic en- 
ergy. This crossing can significantly alter the amplitude 
of the modulation depending on the energy interval cho- 
sen to integrate the differential rate. 

The calculated crossing energies, Ec, (dotted and 
dashed lines) are shown in figure 1, for the case of a 
heavy Dirac neutrino, as a function of the WIMP mass. 

TABLE I. List of parameters used. 

phaio = 0.3 GeV cm"^ 
Vq = 230 km s'^ 
Vesc ~ 570 km s~^ 
cr„ = 270 km s"^ 



Theoretical signal-to-noise ratio 

To extract a signal from the data it is necessary to 
establish the maximum of the theoretical signal-to-noise 
ratio. The total signal as a function of time can be ex- 
pressed in the form, 

S{t) = j ^[V£{t)]dE' ^So + Srn COs(LOi) + 0(^^), 

(2.3) 

where we define So as the unmodulated part, Sm as the 
amplitude of the modulation, uj = 2tt/t {t = 365 d), 
and t is measured from June 2"*^. From the differential 
rates evaluated at June, Sj = S{t ~ 0), and December, 
Sd = S{t = r/2) expressions for Sm and So can be 
obtained. 



So = I [Sj + Sd] ■ (2.4) 



Assuming that Sm ^ 5*0, the theoretical signal-to- noise 
ratio is defined to be. 
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where the product MT is the exposure of the detector, 
namely its mass times the overall exposition time to the 
WIMP flux. This is the function we need to maximize. 



Notice that in (2.5) it is assumed that no spurious back- 
ground is present. A more realistic case will be considered 
in the next section. 



The limits of integration of equation (2.3) are ar bitra ry 
and therefore the theoretical signal-to-noise ratio ( ^.5[ ) is 
effectively energy dependent, 
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In the case of germanium detectors the crossing en- 
ergies, Ec, usually lie between the energy threshold of 
the detector, Ei, and an upper limit, E^ represen- 
tative of the highest energy depositions produced by 
WIMP's. The maximum signal-to-noise ratio, for low- 
mass WIMP's, are obtained when the differential rates of 
June and December are integrated between E > Ec and 
Ef^ since the spectra there differ over a larger fraction of 
the total energy window, {Ei, Ej). For heavier WIMP's, 
however, with larger values of E^ the maximum signal- 
to-noise ratios will be attained when the differential rates 
is integrated in the low energy region, from Ei to E < Ec- 

The theoretical signal-to-noise ratio for 1 kg yr expo- 
sure as a function of the WIMP mass is shown in Fig. 
2; only integrals from E to Ef are plotted. For each 
integration interval the corresponding ratio is calculated 
and plotted at E, while E is varied between and Ef. 
From the figure one can extract, for each WIMP mass, 
the value of E that maximizes the signal-to-noise ratio. 
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FIG. 2. Theoretical signal-to-noise ratios in Ge for 1 kg yr 
for several WIMP masses. The ratios were calculated for all 
energy intervals running from _E to _E/ = 50 keV and plotted 
at E. 



*An upper limit of 50 keV seems adequate since the energy 
deposited in a germanium detector by a 10 TeV WIMP with 
a mean velocity of 340 km s~^ is of that order. 



In a similar fashion, we define E'^^^ as the value of E 
that maximizes {s/n)th when E varies from Ei up to E. 
In Fig. 3 the values of Eg/n and E'^^^ are plotted as a 
function of the WIMP mass. 
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FIG. 3. Values of Eg^„ (^s/n) Gle maximizing the theo- 
retical signal-to-noise ratios when the differential rate is inte- 
grated between E and Ef {Ei and E). 

For the next section a quantitative measure of the an- 
nual effect is needed. We calculate therefore, the percent- 
age ratio of the modulated over the unmodulated signal, 
a = Sm/So, and plot it as a function of the mass of the 
WIMP in Fig. 4. The solid-line is obtained using the 
energy intervals that maximize {s/n)th{E) with relative 
intensities between 2 - 8% (solid line) . In the same figure 
a similar ratio but integrated between an arbitrary min- 
imum energy, Ei = 2.5 keV, and Ef — 50 keV (dashed 
line) is also shown to illustrate the importance of the 
choice of the energy interval. This second curve shows a 
dramatic decrease in the modulation percentage for the 
mass that symmetrizes, with respect to E^ the interval 



WIMP Mass (GeV) 

FIG. 4. Percentage of predicted annual modulation in Ge 
in two illustrative cases. An optimized case (solid line) us- 
ing the energy intervals that maximize {s/n)th{E) and a 
non-optimized one (dash line) where Ei has been arbitrarily 
selected to be 2.5 keV. 
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of integration (« 300 GeV). Table || lists the values of 
Eg/n and the second with three different threshold 

energies, and the corresponding values of Sm/So- 

III. DETECTION SENSITIVITY 

Generalizing the concepts introduced in the previous 
section the total rate under experimental conditions can 
be expressed as, 

Stot ^B + So + 5™ cos(wi), (3.1) 

where B is the rate of spurious background. The experi- 
mental signal-to-noise ratio {Sm ^ So) is then, 

(s/nU, ^ ^VmT « ^^VmT, (3.2) 

where a — Sm/So is as in Table ||. Relation ( |3.2| ) is 
always valid; however, determining the real magnitude 

TABLE II. Energies that maximize the signal-to-noise ra- 
tios in Ge. Depending on the WIMP mass the signal-to-noise 
ratio was obtained integrating the differential rates of June 
and December, dR, in two cases: when dR was con- 
sidered, -Es/n coresponds to the value of E maximizing the 
expected signal-to-noise ratio, and when dR was com- 
puted, -B^/n coresponds to the value of E maximizing the 
signal-to-noise. Ef = 50 keV and Ei = 0, 2, and 3 keV were 
used. 
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of the background implies modeling very accurately sev- 
eral possible sources of spurious pulses and their inten- 
sities (which is rather difficult considering the ambigu- 
ities involved). At this point, one generally assumes 
that Sm is simply a given fraction, a', of the total rate 
B + So and not just a fraction of So', if this is done then 
s/n fa a'^y{B + So)MT. Notice, however, that an ex- 
periment with a large background could drive s/n to ar- 
tificially large values leading to a wrong estimate of the 
sensitivity of the experiment. 

To solve this apparent inconsistency we need a way to 
distinguish between B and 6*0. We propose to use the 
experiment that provides the best exclusion plot as an 
upper bound to the unmodulated dark matter signal, 6*0. 
In other words, we consider the rate of such experiment 
as free of background {B — 0). We call this rate Sl and, 
consequently, for this experiment So = <S'l[]. An arbitrary 
experiment with a total rate higher than that of the best 
exclusion, will then have Stot = B + S^, where B now 
represents the difference in total rates between the two. 
The necessary exposure for this arbitrary experiment, to 
achi eve a level of sensitivity a, can be estimated from 
dJ). Thus, 

- (^)' i ^) ■ 

The best scenario, just mentioned, is that where 7 = 
Stot/ Sl ~ 1, -B ~ 0, and the minimum required exposure 
is therefore. 

Conversely, when Stot/SL ^ 1, longer exposures are re- 
quired to achieve a similar sensitivity. In this case the 
necessary exposure is, 

MT = -/MTm^n ~ ^ MTm^n (7>1)- (3.5) 

To obtain Sl we resort to the exclusion plot, alm^), 
using the relation, 

Je,^„ dE 



'Of course, we do not know for sure if Sl is purely a dark- 
matter signal or still has some spurious background hidden in 
it. The best we can do is to assume that it is pure dark matter 
and estimate the exposure required to detect a modulated sig- 
nal of value a. If no modulation is found, it is either because 
there is none, or because Sl So, and consequently our orig- 
inal assumption of a negligible background was wrong. The 
recipe then calls for a further reduction of the background 
and an increase in the exposure. 
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The differential rate, dR{a)/dE, on the left-hand side 
is calculated as in section II but using the experimen- 
tal exclusion, a{m^)^ instead of the theoretical zero- 
momentum cross section. 



da 

It 



T 

' rr 



F\k). 



(3.7) 



To maintain self-consistency we integrate the differential 
rate using the energy limits and the value of a derived 
in the previous section for each WIMP mass. Thus, the 
predicted signal is renormalized to make it compatible 
with the values of <7{mx) obtained from the exclusion 
plot. 



Solving from equations (B.4) and (3.6), the minimum 



exposure can be expressed as 




(3.8) 
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FIG. 5. Predicted minimum exposures in Ge as a function 
of WIMP mass for three energy thresholds, Ei — 0, 2, and 
3 keV. The curves were derived using the exclusion plot of 
Ref. i. 



In Fig. 5 and Table HI, we show the minimum expo- 
sures (in kg yr) for the lowest rate experiment as a func- 
tion of WIMP mass based on the modulation percentages 
calculated in Table As a reasonable criterion to dis- 
tinguish the signal, we have used that s has to be at 
least 2a larger than the statistical uncertainty or, what 
is the same, that s/n = 2. The figure and table were ob- 
tained using for a{m^) the data from Ref. properly 
rescaled to germanium nuclei. As mentioned in the pre- 
vious section, for low-mass WIMP's, the energy intervals 
maximizing the signal-to-noise ratios tend to be located 
above the crossing energies. For these masses the min- 



imum exposures are obtained integrating relation (3.S) 
between E's/n and Ef. Heavier WIMP's maximize their 
s/n at energy intervals below their crossings. Therefore, 
the three curves corresponding to different energy thresh- 
olds were obtained integrating (|3.8[) between Ei — 0, 2, 



and 3 keV and Notice that the sensitivity decreases 

roughly as a~^. 

For an arbitrary experiment with a total rate Stp t th e 
required exposure can be estimated from equation ( [3.5| ), 
in units of the exposure of the lowest rate experiment. 
The difference between the two can be attributed to a 
spurious background of magnitude B = (7 — 1)Sl- 

An analogous calculation can be done for a Nal de- 
tector using the measured relative efficiencies (quenching 
factors, q{I) = 0.09, q{Na) = 0.30 0) and generalaz- 
ing the formulae of section II to the case of a two-nuclei 
target. Similar results to those of previous sections are 
found with the significant distinction that because of the 
rather low quenching factor of iodine, the values of Ec, 
Es/n, and E'^^^ occur at low energies (< 10 keV). For 
Nal a lower value of Ef was used (30 keV). Furthermore, 
since the iodine nucleus cross section is larger than that 
of sodium, iodine is the main contributor to the total rate 
at low energies (at high energies its nuclear form factor 
strongly diminishes the interaction rate). 

The minimum exposure for Nal as a function of WIMP 
mass is shown in Fig. 6. Due to the large cross section 
of iodine smaller exposures than in the case of germa- 
nium are expected. Nevertheles, as can be appreciated 
from the figure, MTmm is now more dependent on the 
threshold energy of the detector. 
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FIG. 6. Predicted minimum exposures in Nal as a function 
of WIMP mass for three energy thresholds, Ei = 0, 1, and 
2 keV according to the exclusion plot of Ref. [El . 



IV. CONCLUSIONS 

We have performed a quantitative prediction of annual 
modulation rates for Ge and Nal detectors. To this aim, 
we found first the adequate energy intervals where to 
integrate the rate so that the theoretical signal-to-noise 
ratio is at its maximum. Using these optimized theoreti- 
cal rates and assuming that the lowest experimental rate 
achieved so far contains no spurious events we then esti- 
mated the minimum exposure required for the predicted 
modulated rates as a function of WIMP mass. From the 
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results obtained we can draw a few conclusions. To avoid 
artificially large signal-to-noise ratios, an arbitrary exper- 
iment with a total rate, Stot, has to be compared with 
the best experiment, Sl, and should attribute the differ- 
ence in counting rate to the presence of a background of 
the order of i? = (7 — 1)Sl, 7 ~ Sioi/Sl- This arbitrary 
experiment requires an exposure of MT — ^MTmm to 
reach the same level of sensitivity as the best one. Cur- 
rent exclusion plots correspond to exposures in the range 
MT ~ 10 - 60 kg yr for masses > 100 GeV (Ge) and 
between MT ~ 5 — 90 kg yr for masses > 300 GeV 
(Nal) . A strong dependence in Nal of the minimum ex- 
posure with threshold energy is also predicted. In Ge, 
the sensitivity is independent from the threshold energy 
for low WIMP masses. 

A recent claim by the DAMA collaboration jl5| de- 
serves a comment. The best exclusion plot used through- 
out this work corresponds to an experimental mean rate 
of at least, {Sl) « 2 c/kg d (energy region 2-20 keV), ob- 
tained by the same group after rejecting events by pulse 
shape analysis. For the annual modulation analysis, how- 
ever, no rejection is allowed and therefore the correspond- 
ing experimental rate is higher, Stot ~ 1 c/keV kg d x 18 

TABLE III. Rates normalized to be compatible with the best exclusion plot of Bernabei et al. Q and corresponding exposures 
for several WIMP masses in Ge for the intervals of integration of Table II. 



keV « 18 c/kg d. This implies the presence of a spuri- 
ous background that increases the necessary exposure re- 
quired to detect a modulation by a factor 7 = Sioi/Sl ~ 
9. According to our calculations then, the required ex- 
posure 65 kg yr for a 60 GeV WIMP) was not 
achieved by the DAMA group since, at the moment of 
their announcement, they had analyzed a statistics of 
only 12.5 kg yr. It is interesting to note however, that 
our results predict that the DAMA collaboration with its 
115.5 kg detector should be sensitive to low 90 GeV) 
as well as high-mass 1 TeV) WIMPS in approximately 
three years of data acquisition. 

Future improvements in background reduction will cer- 
tainly push further down the limits of dark matter signal 
and at the same time will make the small annual fluctua- 
tions harder to disti ngu ish. Nevertheless, in view of Figs. 
5, 6, and equation (3^), higher-background experiments 
with large masses will still be suitable to explore the low 
and high mass regions. 
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